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Abstract 

We analytically evaluate in the leading logarithm approximation the differen- 
tial cross section for e + e~ — > e + u e u d. We compare our order a A oP s leading-log 
result to the order a 4 a® exact result obtained from the GRC4F Monte Carlo 
program. Finally we use the Gliick, Reya, Schienbien distribution of partons 
in a virtual photon, which incorporates both evolution and nonperturbative 
strong interaction contributions, to obtain better estimates of the differential 
cross section. 
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I. INTRODUCTION 



Events at a high energy e + e~ collider in which only hadrons are seen in the final state and 
the net momemtum of the hadrons transverse to the beam axis is large could be signatures 
for physics beyond the standard model However, these events could also arise from the 
simple standard model process depicted in Fig. |. Here the positron emits a slightly virtual 
photon and escapes down the beam pipe. The hadrons recoil against a high neutrino. 

The cross section for background processes that arise from graphs like that shown in 
Fig. |l| can be computed exactly to order a 4 a° by using the Monte Carlo code GRC4F . 
However, it is useful to have an approximate analytical calculation of the cross section as 
a check, if the calculation is simple enough to illuminate the basic physics. In this paper, 
we provide such a calculation based on the leading logarithm approximation in which the 
photon carries small transverse momentum pt while the quark carries much larger transverse 
momentum k T : ml <C p\ -C k\ <C M^. We compare our order a 4 a° leading-log result with 
the order a 4 a° exact result. 

We also go beyond the a 4 a® leading-log result by incorporating two strong interaction 
effects: the DGLAP evolution of parton distribution functions in a virtual photon giving rise 
to contributions of order J2N( a s In Q 2 /\p 2 \) N , and the phenomenological, nonperturbative 
splitting of the photon into a quark plus anything. To that end, we use the Gluck, Reya 
and Schienbein (GRS) parton distributions in a virtual photon ||. We investigate whether 
these two strong interaction effects are big enough to substantially affect the cross section. 




FIG. 1. Charged current deep inelastic scattering e + e 



e + u e ud. 
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II. FRAMEWORK 



In this section, we set the framework for the approximations we will make and the 
momentum cuts we will use. We first present notation that will prove useful later on in the 
calculation. We will use light-cone coordinates |4j to describe momentum four- vectors 



This set of coordinates is very convenient for symmetric collisions because one of the incom- 
ing particles has almost exclusively plus momentum and negligible minus and transverse 
momentum, whereas the other has almost exclusively minus momentum. 

There are constraints on the momenta of the photon and the quark coupling to the W 
boson. As the transverse momentum squared p\ of the photon falls below the mass squared 
r/ij of the positron, the matrix element squared falls off sharply. We can thus take w? e to be 
an effective lower limit of integration for p\. Its upper limit is the squared veto momentum 
of positron detection, M 2 . This is an experimental limit and will depend on the exact 
configuration of the detector. If the scattered positron has an angle 9 greater than some 
limit #0) it w iH be detected. Thus we demand that \9\ < 6q. In our approximate calculation, 
this amounts to requiring that px < M = 9 0y /s/2. 

A reasonable lower cut-off on the quark transverse momentum squared k\ is the greater 
of m 2 , the mass squared of the p meson (or indeed some other meson such as the n), and 
Pt, i-6. k\ > ma,x(m 2 ,p^). For k\ < m 2 p , the distribution function of quarks in a photon 
is nonperturbative. In the first part of this investigation, we have set this nonperturbative 
contribution to zero. The approximate upper limit on k\ is the virtuality of the W boson, 
Q 2 , which is fixed by the transverse momentum of the neutrino. A plot of these constraints, 
Fig. [|, makes it easy to visualize the phase-space under consideration. 

The approximation we make is to calculate the cross section by including the region in 
the middle of the polygon in Fig. [|, without worrying about what happens at its edges. This 
is reasonable because most events will lie in this region, 



This approximation is called the leading-log approximation, for reasons that will become 
apparent later on. We can split our calculation into three steps, in which the hardness 
increases at each step. 

The first step is to assume that the transverse momentum pt of the photon is much 
smaller than the transverse momentum kx of the quark that couples the photon to the 
W boson (Graph A, Fig. |]). We can then calculate the piece of Fig. [I] dealing with the 
emission of the photon from the positron. The second step is to assume that the transverse 
momentum k? of the quark is much smaller than the transverse momentum of the W 
boson (Graph B, Fig. We then calculate the piece of Fig. |I] in which the photon splits 
into a quark-antiquark pair. The third and final step involves no further approximations 
and is a direct calculation of the cross section for charged current deep inelastic scattering 




Pt = (p\p 2 )- 



(1) 



m l < Pt < k T < Q 2 - 



(2) 



(Graph C, Fig. §). 

The general form of the cross section for the process in Fig. [I] is the following: 
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In k\ 




FIG. 2. Momentum cuts. 



1 2 d 2 P^dP' + d%dl l+ d 2 k' T dk' + d 2 k'^dk" + 

a ~ 2P°2/°2' A ' (2tt) 3 2P'+ (2tt) 3 2/'+ (27r) 3 2fc'+ (2tt) 3 2P+ 

x (2tt) 4 5 4 (P^ + F- P'" - - k'" - jfe""). (3) 

We are now ready to calculate the total matrix element squared. 



III. GRAPH A: CALCULATION OF P y/e 



As just discussed, we will calculate the total matrix element squared |A^a| 2 in three 
successive refinements. First we concentrate on the part of the diagram involving the photon 
(Graph A, Fig. §), which will be the starting point. This is the first step in our three-step 
calculation, as we go from the least hard to the most hard part of the diagram in Fig. [I]. 
|^Vf a | 2 w iH contain a factor |7Wb| 2 that represents the next hardest contribution, while |tV^b 1 2 
will contain the very hardest factor lA^cl 2 - 

We calculate the probability P 7 / e of finding a photon in a positron by evaluating Graph A 
(Fig. D. Avera ging over the initial spins of the positron, we find the total matrix element 
squared to be 

\M A \ 2 = £ Tr {f <p, j) f tfp, i)} 5 ij |^ B | 2 . (4) 

Here we neglect the positron mass compared to px- |A^^| 2 is averaged over the photon 
polarizations and contains the next hardest part which we calculate later on. e(p,i) is the 
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FIG. 3. Graph A. 

polarization vector in the physical A + = gauge of the photon with momentum p^ and 
polarization index %. 

We work in the cm. frame of the electron and positron, with the positron carrying 
only plus momentum and the electron carrying only minus momentum. In this frame the 
momenta in our graphs have components v u = (v + , v~, vt) given by 

P» = (P+,0,5r), 

r = (o,r,o T ), 

Now, e M (p, i) is 

e + (p,i) = 0, e~{p,i) =p T -e T /p + , e 3 (p,i)=5j. (6) 
Using J2i e2 (p, i) — ~ 2, we have 



| ^ A|2= 4(2,) a(1 -fl |l ± (^) !n ^ B|2 (?) 
We recognize the standard splitting function P 7 / e , 



PllM) = I±<Iz«I (8) 

Define p^ ~ p M = (£-P + , 0, 0y). Since the process is very hard, i.e. pt <C kx, the minus 
components and transverse components of p^ are small compared to the plus components. 
In this approximation, 
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FIG. 4. Graph B. 



(2tt) V(P^ + P - P'" - l ,fM - k'" - k"") « (2tt) 4 5 4 (^ + P - - fc'" - (9) 
Then the cross section can be reexpressed to depend on pt and £, 

_2ad 2 p T d£ 1 2 d%dl' + d 2 k' T dk' + d 2 k'±dk" + 

°~ (2n) 2 p 2 T 7/e ^2p°2Z°2' 1 (2tt) 3 2Z'+ (2vr) 3 2A;'+ (27r) 3 2fc"+ 

x (2vr) V(^ + P - /' M - ^ - fc""). (10) 



IV. GRAPH B: CALCULATION OF P g/7 

We now calculate the probability of finding a quark in a photon. This will give us |7Wb| 2 , 
in terms of the final hard matrix element squared |A^c| 2 - There is also a graph B' in which 
the quark is replaced by an antiquark. In this case, everything is the same except the hard 
matrix element |A^c| 2 is replaced by a different matrix element |A^c'| 2 - 

Summing over colors and averaging over spins, the matrix element squared for Graph B 
(Fig. |) is 

2 2 

I-MbI 2 = N c^y\ E Tr {H ¥ rfp, i) Y' rffr i) V} ■ (11) 

Again, the mass of the quark is negligible compared to kf. The matrix Ti contains the final 
contribution that we will evaluate using Graph C (Fig. [5]) later on. For now, we leave it as 
it is. Let Jj —ft ^(p,i) ft" /{p,i) ft- Since the partons involved are going much faster in the 
plus direction compared to all other directions, we can make the approximation Ju ~ L + 7~. 
Then, 
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I^B| 2 = iV c ^I^Tr{H 7 -}L^. (12) 
Let k^ = (k + , 0, T ). We use it to rewrite the above, 

L+ = Itt { 7 +l} = Ity { 7 + K rfp, i) K'7(p, i) K} ■ (13) 

The hard scattering matrix element squared, which we leave unevaluated for now, is |7Wc | 2 = 
|Tr 1 7i J^j, and the matrix element squared is in the simple form 

2 2 

I-MbI 2 = N c ^- 2 ^ ~Tr { 7 + K rfp, i) r rfp, i) K} |M;| 2 . (14) 

In the same cm. frame as earlier, but now with the approximation p M — > p M described 
above, 

p" = (p + ,0,0 T ), 

•"'(('-^^-^ (15) 
After doing the trace algebra, the matrix element squared is 

I-MbI 2 = 4(2 J )a(1 -^ iVc (, 2 + (i - ;?} \M C \\ (16) 

We recognize the standard splitting function P 9 / 7 , 

P ff/7 (*)=tf C {^ + (l-*) 2 }. (17) 



In the same approximation as before, the minus and transverse components being small 
compared to the plus components, the delta function in the cross section can be expressed 
as 

(2tt) 4 (5 4 (^ + P - Z'" - k >tx - « (2tt) 4 (5 4 (^ + F - - k'^). (18) 

Let x = k + /P + , which implies k + = (x/(,)p + and z = x/£. Then, using the result we 
obtained for |.Mb| 2 , the cross section can be written as 



1 . ,„ d 2 lUl' + d 2 kUk' + 



X 



^c| 2 ^ L,„ + , ^ , /+ (2tt)V(P + P - Z'" - (19) 



2fc°2Z°2' ' (2tt) 3 2/'+ (2tt) 3 2A;'+ 



which makes it obvious that the cross section is the convolution of a hard-scattering cross 
section with a piece that we identify as the distribution function of a quark in a positron. 

A brief description of the required distribution functions is relevant. Graph A (Fig. [3]) 
involves the distribution function of the photon in the positron, and Graph B (Fig. |4]) the 
distribution function of a quark in a photon. As we combine the two graphs, we will require 
the distribution function of a quark in a positron. We will perform the trivial angular 
integral to write d 2 px = ^dp\. We now define the distribution functions described, 

f M -— r'Mp ( \ 



2-K J ml Pt 

f qh {x,Q ,p T ) = -~r , a 2 , -Tr p qft( x h 

f,/^Q 2 ) = ^f%MP,/^mJ M ' d 4l Q ' l2 (20) 

Jx c; Jm% Pj> Jmaj[(py,mj) tvrp 

We first evaluate the convolution of the splitting functions, 

t f *7/.(0JVr(*/0 = ^ ((! " *) + 2 (! + ^) Inx + • (21) 

Next, we compute the momentum integral, either by inspection of Fig. |2|, or by explicit 
calculation, 

M 2 dp 2 Tf Q 2 dk 2 T Q 2 m 2 Q 2 M 2 1 2 M 2 



I ml p\ imax(p|, m 2) k\ Til 2 Till M 2 TTl 2 2 Wb 2 

Therefore, the distribution function of a quark in a positron (or, equivalently, in an electron) 



is 



a 2 e 2 ( 4 (1 — r 3 l 

f q/e (x, Q 2 ) = — ^N c (1 - x) + 2(1 + x) Inx + -- 



(2tt) 2 ° V 3x7 

x (\n ^- In ^5 + In ^— In — + - In 2 —) (23) 
\ m 2 ml M 2 m 2 2 m 2 J 

If we use phenomenological distribution functions / g / 7 of quarks in virtual photons (e.g. 
GRS distribution functions), we can evaluate f q / e as 

f q /e(x,Q 2 ) = ^z\ Xd J ^ ^fP,/e(x/Of q /^Q 2 ,p 2 T)- (24) 
i/T 15 Jmf Pj> 

We will use this convolution in a following section when using the GRS distribution functions. 
We label the hard scattering cross section da, 

1 r\ 2 V r\ 2 V rlh' + 

da = -J—\Mc\ 2 t \J % A 2itY5\W + l* - h'»). (25) 
2fc°2Z°2 (2tt) 3 2/'+ (27r) 3 2fc /+ v J K ' K ' 

The cross section da is a convolution of the distribution function with the hard scattering 
cross section 



da 



J dxf q/e (x)da. (26) 
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V. GRAPH C: CHARGED CURRENT DEEP INELASTIC SCATTERING 



The third and final part of the calculation is to compute the charged current deep inelastic 
scattering cross section, Graph C (Fig. ||). This final calculation is quite standard and 
involves no further approximations. It is a direct calculation of the Feynman diagrams in 
Fig. It involves the W boson coupling constant gw = V&ra/ sin 9w and the element of 
the quark mixing matrix V u d that relates the up quark to the down quark. 




Graph C Graph C 

FIG. 5. Charged current deep inelastic scattering. 
Averaging over initial spins, the matrix element squared for Graph C (Fig. |||) is 



_ 2s 2 (2n) 2 a 2 \V ud \ 2 



where s is the cm. energy squared of the electron-quark system. 

The calculation of the matrix element squared for an incoming antiquark instead of a 
quark is nearly identical, and is depicted in Graph C, Fig. The matrix element squared 
is 

lMc ' 1 "^MW ( H) ' {) 

where y = q ■ P/l ■ P. 

We are now ready to evaluate the hard scattering cross section. We first rewrite it in a 
more suitable form, 

1 rP-V rll' + 

dcr = d % %-5{k n )\Mc\ 2 . (29) 

2A;°2/°2 (2vr) 2 21'+ y n C| y ' 

Now, 

k' 2 = (k + qf = 2xq ■ P - Q 2 = Q 2 (x/x bj - 1), (30) 

where = Q 2 /2q ■ P. We express the momentum /' of the outgoing lepton in terms of its 
transverse momentum V T and its pseudo-rapidity 77', 
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Z'" = (l' T e*/V2, l' T e-*/V2, l' T ) , (31) 

so cZ/ /+ = l' + di]'. The azimuthal angle integral is trivial, rf 2 /^ = 7r<i/^. In terms of the cm. 
energy of the hard process, k°l° — s/4. 
The hard cross section is then 

d * = T^-^ dl ?dv'S(x - x bj )\M c \ 2 . (32) 

107TS L> 

Using the squared matrix elements |.Mc| 2 and |.Mc'| 2 we calculated previously, and the fact 
that Q 2 = ys, we obtain 



~ vra \V ud \ x ,,/2, it, 

da e - u ^ Ved = , 4 -dl T dr] 5{x - 

^-^• fl = 2sin^(Q^+^)2 y d ^ (g - - ^ (33) 

VI. DIFFERENTIAL CROSS SECTIONS 

We have computed all the pieces necessary for the calculation of the cross section, 

da = J dx {f u /e(x)da e - u ^ Ued + fd/ e (x)da e - d ^ Uea } . (34) 



x bj ) 



The result is 



dl' 2 dr)> %<k 6 W {Q 2 + M 2 ,) 2 \2 m 2 m 2 m\ M 2 m 2 ) 

x [x{\ -x) + 2x(l + x) lnx + ^(1 - a: 3 )) (^±3^^ . (35) 

It is often more convenient to rexpress this in terms of l' T rather than l®, 

d ° (36) 



dl' T drf v "dlfdrf 

The final step is to change variables from x, y, Q 2 to rj', l' T , s where s is the cm. energy 
squared of the e + e~ pair, not to be confused with the hard scattering cm. energy squared 
s. In the laboratory frame, neglecting m e in the change of variables, 



p"= (7^2,0, o r ), 



M^^?4 (37) 
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The change of variables is 



qp = -{l-l') 2 = e lf t r ^, 
Q 2 e*l' T 



x 



2q-P y/s-e-Jlir' 



q-P ILe-* , s 

v = TT = l - s vr- (38) 

The physical region is where x < 1: 

2^ cosh r/ < yfs. (39) 

We now have a simple approximation to the cross section that can be compared with 
the output of Monte Carlo programs. We check our cross section against GRC4F 0. Later 
we introduce a refinement not included in the purely perturbative Monte Carlo calculation. 
From GRC4F, we obtainQ the momentum four-vectors of all final state particles for 3666 
events of the type we are interested in corresponding to a luminosity L = 100 pb~ . The cm. 
energy is \fs = 189 GeV and the veto momentum squared was taken to be M 2 = 5 GeV 2 . 
To be consistent with the choice made in getting the GRC4F results, we pick the cut-off 
for the pertubative treatment of our parton distributions to be twice the pion mass, i.e. 
kj< > max((0.28 GeV) 2 ,PT). We make a histogram of cross section da/dl' T versus transverse 
momentum of the neutrino l' T . Next, we compare this GRC4F result to the approximate 
cross section obtained by numerically integrating the differential cross section we calculated 

da fi» da 

d% = L v ^ dr] di^/- m 

The range of pseudo-rapidity follows directly from the physicality condition in Eq. 



v - = log 2v7i • {) 

We plot this theoretical curve over the histogram in Fig. || Taking into consideration the 
multiple approximations we made, we see that our calculated cross section agrees quite well 
with the Monte Carlo GRC4F, except at low l' T , where the leading logarithm approximations 
used in the approximate calculation do not apply. 

We follow a similar procedure to obtain a graph of da/drj'. We bin the events in pseu- 
dorapidity rj' of the neutrino, from rj' = —10 to rj = 10, which is a reasonably large enough 
range to allow the cross section to fall off sharply at its edges, and we plot a histogram. 
Finally, we numerically integrate the differential cross section we calculated, 

da /■^Tmax ,„ da 

— = / dl' T — . (42) 

di Jo T dv T di v ; 



We thank Alain Bellerive of the OPAL collaboration at CERN for providing us with the GRC4F 
Monte Carlo events. 
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da/dl' T , pb/GeV 
1.2 




Z^, GeV 



FIG. 6. Differential cross section da/dl' T versus transverse momentum of the neutrino. Our 
theoretical line is superimposed on results from Monte Carlo simulations using GRC4F with 3666 
events. 



From the boundary conditions of the physical region, the integration range of l' T is from 
to l' n 



Tmax' 



l'n 



Tmax 



2 cosh 7]' 



(43) 



The histogram and the theoretical curve are plotted in Fig. |7|. Once again, we see 
that our approximate calculation and GRC4F agree except for large negative rj', where the 
approximations we used break down. 

As a final check, we compute the total cross section a. The GRC4F cross section is just 
the total number of events of this process divided by the luminosity, so 0"grc4F = 36.7 pb. 
We calculate the total cross section from our result by either integrating da/dl' T or da/drf 
over their entire ranges < l' T < y/s/2 and — oo < rj' < oo. The cross section we obtain by 
either method is cr t h = 30.7 pb. This is in reasonable agreement with the GRC4F result. 

In conclusion, we see that the cross section we calculated analytically from first principles 
in the leading log appromation is not far off from the GRC4F cross section in the region of 
large l' T and not too negative rj', where the momentum transfer from the electron is large. 



VII. REFINEMENT USING THE GRS PARTON DISTRIBUTIONS 

So far, our calculation has been at order a°, using pointlike parton distribution functions 
of quarks in a photon, without any evolution. The GRC4F result was exact at order a 4 a°, 
and our approximate calculation was a leading logarithm calculation at order a 4 a°. In this 
section, we extend the calculation to include strong interaction effects. We continue to work 
in the approximation that virtualities are strongly ordered (the leading log approximation) 
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da/drf, pb 




-6 -4 -2 2 4 



FIG. 7. Differential cross section da/drj' versus pseudo-rapidity of the neutrino. Our theoretical 
line is superimposed on results from Monte Carlo simulations using GRC4F with 3666 events. 

since we have seen that this approximation works well for large and rj' ~ —3. We inves- 
tigate two strong interaction effects. First, the emission of gluons from the quark before 
it is struck can lead to contributions of order J^n^s^Q 2 /\p 2 \) N , that tend to lower the 
cross section. We can account for these potentially important contributions by using parton 
distribution functions / g / 7 that have evolved according to the DGLAP evolution equation. 
Second, for \p 2 \ $ m 2 , there can be important nonperturbative hadronic contributions to 
fq/^. This "vector meson dominance" component can be modelled by using the parton 
distribution functions in pions 0, which are determined empirically. We can account for 
both these effects by using the recently published || parton distribution function for virtual 
photons by Gliick, Reya, and Schienbein (GRS). 

We define the virtuality of the photon to be — p 2 . In the calculation that follows, we will 
make the approximation p 2 = —p\. The GRS parton distribution of a virtual photon has 
the form 

f$* (x, Q 2 , -P 2 ) = ff h (x, Q 2 , -P 2 ) + f q ^(x, Q 2 , -p 2 ) 
= ff h M 2 ^ v 2 ) 

+ (1-^2)2 (<%M*>&) + ^\{G 2 U - G 2 d )f s/7T (x,Q 2 )) , (44) 

with G 2 U ~ 0.836, G 2 d ~ 0.250, and 5 U = —1, 5 d = +1. The two pieces of this parton 
distribution are the pointlike piece /^J 7 and the hadronic piece fqf^- The latter is a function 
of the parton distribution in a pion f q / n and of the strange sea distribution f s / n , weighted 
by factors of order unity and a overall factor (1 — p 2 /m 2 )~ 2 which turns on sharply as the 
virtuality — p 2 of the photon goes below the mass squared, m 2 , of the p meson. Here ir 
refers to the neutral pion tt°. The nonperturbative hadronic part is seen to start playing 
a role in the overall parton distribution at low photon virtuality, where the photon can be 
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da/dl' T , pb/GeV 
1.2 




— - — - total 
a 4 a° leading-log 

— — — - pointlike component 
hadronic component 



l' T , GeV 



FIG. 8. Differential cross section da/dl' T versus transverse momentum of the neutrino. Our 
theoretical a 4 a® leading-log calculation and the contributions due to the pointlike component of 
the GRS parton distribution, the hadronic component of the GRS parton distribution and the total 
GRS parton distribution are superimposed on a 4 a^ results from Monte Carlo simulations using 
GRC4F with 3666 events. 



naively thought of as a vector meson. The distributions f^, , f q /„ and f s / n are specified 
parametrically by Gliick, Reya, and Schienbein. We compute the distribution of quarks in 
an electron using the convolution 

/* G /e S (*> Q 2 ) = Tj X T /T ^^(-/O/Sf Q\V%). (45) 

We perform the calculations of the differential cross sections as described earlier, using 
the new GRS distributions. The results are displayed in Figs. || and [| 

We notice a few interesting features of the differential cross sections calculated us- 
ing the GRS distributions. The cross section due to the pointlike GRS distribution has 
the general shape of the cross section we calculated using our unevolved pointlike dis- 
tribution but is somewhat smaller. We should mention that the lower cut-off on kx, 
k\ > max((0.28 GeV) 2 ,p^), we use to match the GRC4F results is lower than the corre- 
sponding cut-off used by GRS, max(0.26 GeV 2 ,PT). If we set our cut-off equal to the GRS 
cut-off, the difference between our unevolved pointlike distribution and the GRS pointlike 
distribution diminishes, but does not vanish. One conjectures that this difference is due to 
the evolution built into the parametrization of the GRS pointlike distribution. The contri- 
bution due to the hadronic piece of the GRS distribution is fairly soft, as one would expect. 
The differential cross sections da/dl' T (Fig. |Sp and da/drj' (Fig. |^) are in reasonably close 
agreement with the GRC4F Monte-Carlo data, except once again at low V T and very negative 

7]'. 

This check of our calculation and of the GRC4F Monte-Carlo leads us to believe that 
the exact result of the differential cross sections should lie reasonably close to the total GRS 
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da/drj', pb 




FIG. 9. Differential cross section da/drj' versus pseudo-rapidity of the neutrino. Our theoretical 
a 4 a® leading-log calculation and the contributions due to the pointlike component of the GRS 
parton distribution, the hadronic component of the GRS parton distribution and the total GRS 
parton distribution are superimposed on a 4 a^ results from Monte Carlo simulations using GRC4F 
with 3666 events. 

curves in Fig. [8] and |9], excluding perhaps the very low l' T and very negative rf regions. 
One estimate of the error would be the difference between the GRS curves and the a°a 4 
leading- log curves. 

We conclude that the GRC4F data is for the most part within the errors of our cal- 
culation. This gives us confidence both in our calculations and the GRC4F Monte Carlo. 
Finally, we have gained some insight into the physical content of the different parts of the 
cross section. 
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